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Abstract. - The strictly nonrelativistic isospectral scheme based on the 
general Riccati solution and Darboux transformation function correspond- 
ing to excited states is presented on the following examples: the harmonic 
oscillator, the square well, and the hydrogen atom. 



PACS ll.30.Pb - Supersymmetry. 



Supersymmetric quantum mechanics (SUSY QM) Q is based on Riccati 
equations of the form 

+y' = Vi(z) + e , (1) 



dz 

where Vi{z) is the initial Schrodinger potential and e < is the so-called 
factorization energy that for convenience we shall absorb in the potential, 
i.e., we shall redefine the initial potential as Vi = Vi -|- e. The SUSY partner 
Riccati equation will be 

-^+y^ = V2(.), (2) 
dz 

where the partner potential V2{z) is Darboux 'isospectral' with respect to 
Vi{z), i.e. 

V2{z)=Vi{z)-2D\\n{^{z))] , (3) 

where D = and t/; is a particular solution of the Schrodinger equation 
D^V ~ ^iV' = 0- The only difference between the spectra of V2 and Vi 
is that the energy level corresponding to the wavefunction ip{z) used as 
transformation function is missing from the spectrum of V2. 
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In 1984, Mielnik introduced a SUSY QM scheme based on the general 
Riccati solution with application to the harmonic oscillator case @. In 
the latter approach, one gets a one-parameter family of Darboux strictly 
isospectral potentials with respect to Vi given by 

Vi iz;X) = Vi (z) - 2D^ [ln(/(z) + A)] , (4) 

where I{z) = Jj^ ip'^{x)dx. Moreover, 



Ji yj^{x)dx + A 

is the modulated Schrodinger mode implied by this scheme in which the Ric- 
cati integration constant A is kept as a free parameter. The inferior limit / 
of the integral I{z) is zero for the radial problems, — oo for the full-line prob- 
lems, and —L/2 for the square well case (see below). The factor ^/X{X + 1) 
is a normalization constant of the strictly isospectral SUSY modes that in- 
terestingly does not depend on the quantum numbers. In general, to get 
continuous solutions ^l'{z,X), one should avoid a zero denominator. This 
leads to conditions on the possible values of A. On the other hand, if one 
works with polynomial solutions, which is the usual case for the discrete 
spectrum of an exactly solvable quantum problem, there will be singulari- 
ties in the logarithmic derivative of ip. SUSY partner potentials based on 
the n'th excited state of a Schrodinger discrete spectrum problem split in 
n + 1 branches separated by the n singularities of the logarithmic derivative. 
Such problems have been considered by Robnik Q. 

In the following, we focus on Robnik's results from the perspective of 
the strictly isospectral scheme, namely instead of working with Eq. (^) as 
Robnik did, we make use of Eqs. and (|5|). This is also equivalent to say- 
ing that we work with factorization energies e = e„ in Eq. (|l|) and ijj = ipn 
as transformation function in Eq. (Q), where e„ is an eigenvalue of the dis- 
crete spectrum of the exactly solvable problem and is the corresponding 
eigenf unction. Our applications are the following. 

(i) The one-dimensional harmonic oscillator 

The initial potential is 

Vi{z)=z^ -{2n + l) , (6) 
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for n G 0,1,2,.... The wave functions are given in terms of the Hermite 
polynomials in the form: 



V'n = , Hn{z) exp(-zV2) (7) 

^2"n!^ 

We present plots of Eqs. (^) and (|5|) corresponding to n = 2 and n = 3 
in Figs. 1^ and 0, respectively. The examination of the plots shows that there 
are no singularities in the one parameter family of potentials for the allowed 
range of the A parameter. Instead, there are n + 1 wells, where n is the 
number of nodes of the excited state on which the strictly isospectral SUSY 
construction is based. Thus, we have at our disposal a definite method of 
producing multiple-well oscillator potentials of exactly solvable type. 



(ii) One-dimensional square well of length L 
The initial potential is: 

I Vi(z) = +00 otherwise , 

where n belongs to the set n = 1, 2, .... The wave functions within the well 
are: _ 

/2 TITTZ 
— cos(— — ) for odd n, (9) 

and: 

[2 mrz 

ipn = y y sin(— — ) for even n . (10) 

Plots of the strictly isospectral formulas for this case are given in Figs. ^ 
and for an odd case and an even case, respectively. 



(iii) The three-dimensional Kepler problem 
The initial potential is 

for rir- from the set = 0,1,2,..., and where the centrifugal barrier has 
been also included. 
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Figure 1: The strictly isospectral potential and the deformed wavefunction 
(multiplied by 10 for the sake of visualization) obtained using the n = 2 
excited state of the harmonic oscillator as Darboux transformation function. 
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Figure 2: The strictly isospectral potential and the deformed wavefunction 
(with the same scale factor of 10) obtained using the n = 3 excited state of 
the harmonic oscillator as Darboux transformation function. 
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Figure 3: The strictly isospectral potential and the deformed wavefunction 
X 10 obtained using the n = 1 excited state of the square well as Darboux 
transformation function. 
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n=2 




Figure 4: The strictly isospectral potential and the deformed wavefunction 
X 10 obtained using the n = 2 excited state of the square well as Darboux 
transformation function. 
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The radial wave functions to be used as Darboux transformation functions 
are 

where z = 2r/n and the Bohr radius has been chosen as unit. 



In Figs, ^-g and |^-^ we present the result of the strictly isospectral 
construction based on the hydrogen radial wavefunctions R^q and re- 
spectively. For the employed range of the A parameter, there is a small 
strictly isospectral effect extending in the Rydberg region till about 20 Bohr 
radii. On the other hand, although not shown in the plots, we can report 
that in the z = nuclear region the efect can be as big as 25%. 

We recall that for this case a similar approach has been used by Fernandez 
However, his procedure cannot be applied to the I = cases, whereas the 
way we use the strictly isospectral Darboux scheme covers the full hydrogen 
spectrum. 

Finally, we notice that it is still an open issue to what definite physical 
situation the strictly isospectral scheme does apply. There are hints in the 
literature indicating that A can be used as a measure of confining effects 
on the quantum spectra. 
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Figure 5: The percentual difference 100(Vi(z; A) — V\{z))jV\{z) between tlie 
strictly isospectral Coulomb potential and the Coulomb potential when the 
hydrogen radial state i?3o is used as Darboux transformation function for 
four values of the parameter A. 
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Figure 6: The difference between tfie probability densities ^/^^(z; A) — ip'^{z) 
for the same case as in Fig. 5. 
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Figure 7: The difference between the strictly isospectral Coulomb potential 
and the Coulomb potential Vi{z; A) — Vi{z) when the hydrogen radial state 
i?2i is used as Darboux transformation function for the same four values of 
the parameter A. 
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